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Abstract. In this review, I focus on the structure and evolution of the coro- 
nal magnetic fields modelled from observations. The development of instruments 
measuring the photospheric and chroniospheric magnetic fields with a high spa- 
tial and time resolutions allows us to improve the modeling of the coronal fields 
based on extrapolation and evolution techniques. In particular, I detail the 
advance modelling of quiet-Sun areas, active regions and full-disc evolution. I 
discuss the structure of coronal magnetic features such as filaments, sigmoids 
and coronal loops as well as their time evolution and instability. The complexity 
of the coronal field and the origin of open flux are also investigated in these dif- 
ferent areas. Finally I discuss the future improvements in terms of instruments 
and models required to understand better the coronal field. 



1. Introduction 

The structure of the solar corona is organised by its magnetic field. The plasma 
/3, the ratio of the plasma pressure to the magnetic pressure, is less than 1 from 
the bottom of the corona to about 2.5 solar radii ( Gary 200 ll ). Unfortunately 



to date, it is not possible to obtain a reliable measurement of the full magnetic 
vector in the corona, whilst it is routinely observed in the photosphere and 
chromosphere. Waiting for improved spectropolarimetric observations based on 
coronal lines, the coronal magnetic field has thus to be derived from physical 
assumptions relying on photospheric measurements. 

The physical assumptions to model the coronal field depend on the areas of 
the Sun considered. Three different parts are often distinguished: (i) the active 
regions, (ii) the quiet Sun, and (iii) the full Sun. The active regions on the 
photosphere are regions of strong magnetic field which can form sunspots. The 
time evolution of active region is 10-15 min by considering an average Alfven 
transit time along a loop of 200 Mm (except during flaring activity) . The spatial 
scale is typically 300-500 Mm. Observations in a broad range of wavelengths 
show coronal structures such as filaments /prominences, sigmoids, loops. Active 
regions are the sources of eruptions such as flares and Coronal Mass Ejections 
(CMEs) due to the large amount of magnetic energy stored. The quiet Sun is 
the area outside active regions (including coronal holes). The typical spatial and 
time scales are imposed by the structure of granules (20 Mm and few minutes). 
The full Sun encompasses both the quiet Sun and active region areas. At the 
current spatial resolution, the time of evolution is imposed by the evolution of 
active regions for instantaneous magnetogram and by the differential rotation 
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for synoptic magnetic maps. The models are imposed by the physical conditions 
of the different regions and constrained by the observations. 

In this review, I only focus on the models developed to determine the dif- 
ferent structures of the coronal magnetic field from photospheric observations 
with a special emphasize on nonlinear force- free models. 



2. Force- free models 

2.1. Magnetic field extrapolations 

Magnetic field extrapolations consist in computing the coronal magnetic field 
assuming an equilibrium state and using the distribution of the magnetic field 
observed in the photosphere or chromosphere as boundary condition. This is a 
static model of the corona. In the corona, three main forces act on the plasma: 
the plasma pressure gradients, the gravitational force and the magnetic forces. 
The equation governing the equilibrium is then: 



- Vp + {V A B) A B + pg = 0. 



(1) 



Several main assumptions are thus defined depending on the time and spa- 
tial scales to describe and, most importantly, on the magnetic field measurements 
available. It is worth noticing at this stage that none of these assumptions can 
describe the real physical nature of the corona as plasma flows have an important 
role. Nevertheless the study of magnetic equilibria remains a key to understand 
better the complexity of the coronal magnetic field and, to date, this is the most 
reliable method to access the 3D coronal field from observations. 

Three main assumptions are currently in use to extrapolate the magnetic 
field into the corona: potential field for which no electric currents (or curl of 
magn etic field) are present in the configuration ( Schmid^ll964l : ISemel &: Ravrole 
19681 ). the linear force- free field in which the electric currents are parallel to the 
magnetic field li ne and the coefficient of proportionnality is the sa r ne everywhere 
in the volume ( Nakagawa Sz Raadul Il972l : IChiu Sz HiltonI Il977l : lAlissandrakis 
19811 : ISemellll988l : lGarvlll989l l . and the nonlinear force-free field in which the co- 
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efficient of proportionnality var ies from one field line to an other (e.g., 
19581 : ISakurail Il98ll : lAlvl 1 19841 ) . The latter assumption is the most realistic 
and most advanced technique in use. Nonlinear force-free extrapolation tech- 
niques can be classified depending on the boundary c ondition s they use or their 
numerical schemes: optimisation (Wheatland et al.l [2000: Wiegelmann 20041: 
Wiegelmann fc Neuki rch'2006': Wiegelmann et al. 2006, 2008; Tade sse et al.ll2009l '). 
Grad Rubin (iGrad &: Rubin.1958.: .Sakuraiiil98L: .Al v 1989; Amari e t al.lll997l. 
1999| :| Wheatlandll2004l:lAmari et al.ll2006l :llnhester fc Wiegelmann 200 6: IWheatlandl 
2006 , 20071: IWheatland fc Regnicr" 2009l') Tevo lutionnary techniques (IMikic McClvmont 



1994 ). magneto -frictional fiYang et al.lll986l: Ivan Ba llegooi ien et al 



cal integration (IWu et al.lll99(]l: iDemoulin et al.lll 992: Son g et al. 
ary integrals (Yan fc Sakurail '2000: Yan fc Li 2006: Valori et al 



20001 ) ■ verti- 
20061) . bound- 
20051 ). In re- 



cent reviews (jSchriiver et al. 2006 : iRegnier. ,2007 : .Wiegelmann 20081 ) . the pros 



and cons of the different numerical schemes are discussed and compared us- 
ing semi- analytical models or observations. To determine nonlinear force- free 
configurations, the boundary conditions are the vertical or radial component of 
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the magnetic field, and either the distribution of a in one polarity or the two 
transverse components of the magnetic field in both polarities. The first set 
of boundary co nditions corresponds to a mathern atically well-posed boundary 
value problem (|Grad fc Rubinlll958l : ISakurailflQSll l. 

Recently, more sophisticated assumptions have been developed to improve 
the physical content of the above models: the magnetohydrostatic model which 
takes into a ccount the plasma pressure grad i ents and/or the gravitational force 
jLow'igSS'; Bogdan fc Lowlll986l : lLowlll99ll : lNeukirchlll99 5': Wicgclmann et al 
2007: Kuan et al..,200§'). an d non force-free models (|Hu DasguDta.200& . ,20o" 
Hu et al.ll2008l : lGarvll2009l ). 



The force-free reconstruction is applied to an observed magnetogram at a 
given time and without a priori on the structure of the coronal field. Several 
other methods have been devel oped to construct force - free equilibria adding 
constraints from observations. In Ivan Ballegooiien ( 20041 ) . a weakly twisted flux 
rope is inserted into a potential field configuration and then relaxed to a non- 
linear force-free state. The flux rope insertion model is constrained by chromo- 
spheric or coronal observations. Unlike the nonlinear force- free reconstructions 
mentioned above, the flux rope insertion model only requires the vertical or ra- 
dial component of the magnetic field measure d on the ph otosphere. In addition, 
a weakly magn e tohyd rostatic model based on lLowl ( 199ll ) has been developed by 
Aulanier et all (|l999l ) imposing an a priori external bipole allowing the existence 



of twisted flux bundle. 



2.2. Magnetic field evolution 

In order to follow the evolution of the solar corona, two approaches can be 
followed. 

First, a time series of equilibria can be constructed from observed magne- 
tograms assuming that the time of evolution of the coronal structures is slow 
enough compared to the reconnection time and the Alfven transit time. The 
method does not consider the history of the region as the magnetograms are 
treated independently. Nevertheless, part of the history of the region is included 
in the electric currents for nonlinear force- free models. The techniqu e of succe- 
sive force- free equil ibria as been applied bv iHevvaerts Sz Priest ( 19841 ) for linear 



force-free fields, bv lRegnier Canfield ( 20061 ) for nonlinear force- free fields. 



Second, the flux transport model is use d to describe the long-term evoluti on 
of the solar corona during a magnetic cycle (iMackay van Ballegooiienll2006al lbl: 



or tne solar corona aurmg a magnetic cvcie (liviacKay &l van .Daiiegooiienllzuubal JDl: 
Yeates et al.ll2007l . l2008al ]bl: lYeates Mackavll2009l '). The flux transport model is 



twofold: (i) evolution of the photospheric magnetic field, (ii) magneto-frictional 
relaxation to a nonlinear force-free equilibrium. The photospheric boundary 
conditions (usually synoptic maps) are evolved in time by including in the in- 
duction equation the effects of differential rotation, meridional fiow and surface 
diffusion. Those three effects have different characteristic times: 0.25 years for 
the differential rotation, 2 years for the meridional fiow and 34 years for the sur- 
face diffusion. Once the photospheric field is determined, the coronal magnetic 
field is derived from the magneto-frictional relaxation method allowing the mag- 
netic configuration to relax to a nonlinear force-free state. This method takes 
into account the history of the magnetic field during a cycle and also includes 
an automatic procedure of emerging magnetic bipoles based on the best match 
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with coincident observations. The flux transport model is used to describe the 
large scale structure of the corona as a nonlinear force-free field, and thus gives 
physical insights different from the large-scale potential field commonly in use. 
As this review focuses on the nonlinear force-free modelling of the solar corona, 
I will omit the magnetohydrodynamic (MHD) models of coronal fields. 



3. Structure of the solar corona 

The above techniques to derive the force-free nature of solar corona have been 
applied to different magnetic regions: (i) the quiet Sun, (ii) the active regions, 
and (iii) the full Sun. 



3.1. Quiet Sun 

The quiet Sun is a misnomer. The evolution of the quiet Sun has a character- 
istic time of the granule evolution of few minutes. Consequently, lots of erup- 
tive events are observed continuously. To model the magnetic field in observed 
quiet-Sun regions, the model used is the potential field because the magnetic 
field measurements are mostly provided by the line-of-sight component of the 
magnetic field. Nevertheless recent development in instrumentation shows the 
possibility to measure reliably the three components of the magnetic field with 
a great accuracy. 

The quiet-Sun magnetic field, the so-called magnetic carpet, has been mod- 
eled as a potential fi eld defining the polarities as point sources. This is the 
point charge method ( Schrijver Sz Title l2002l : ILongcope et al.ll2003l : IClose et al 



l2003l 1. These models are based on SOHO/MDI line-of-sight magnetograms as 
boundary conditions. SOHO/MDI has a mod erate spatial resolut ion of 1.98 
arcsecond and a time cadence of at most 1 min. iReg nier et all f|2008l ) have com- 



puted the potential field of a quiet-Sun region observed by Hinode/SOT/NFI 
with a spatial resolution of 0.16 arcsecond. This model does consider a con- 
tinuous distribution of the magnetic field on the photosphere. In addition, as 
revealed by previous work, the complexity of the qu iet-Sun magne t ic fiel d lies 
near the photospheric surface (below 5 Mm), therefore iRegnier et al. ( 20081 ) have 



implemented a stretch grid along the vertical axis with a very fine grid near the 
bottom boundary in order to resolve the nonlinearities of the magnetic field. 
The authors revealed that the complexity of the magnetic field (defined as the 
number of null points) is concentrated in the photosphere and the chromosphere 
(below 3.5 Mm) whilst the corona above a quiet-Sun region is not complex. 
By measuring the vector magnetic field, Hinode/SOT has successfully measured 
the magnetic field in coronal holes, and thus shows the almost unipolar nature 
of coronal hole: the small bipoles being connected at low height in the chro- 
mosphere or the bottom of the corona. Consequently, the open magnetic fiux 
responsible for the fast solar wind has a strong latitudinal dependence as the 
quiet-Sun magnetic field becomes more and more unipolar from the equator to 
the poles. 

The structure of the coronal field in the quiet Sun is complex and dynamic. 
The magnetic field evolves on the time scale of a granule. 
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3.2. Active regions 

The potential field is a minimum of magnetic energy for a given distribution of 
the vertical or radial magnetic field component on the photosphere. Therefore 
there is no free magnetic energy, no shear and/or twisted field lines in a poten- 
tial field configuration. For these reasons, the nonlinear force-free field is more 
adequat to describe better the nature of the corona as it contains free magnetic 
energy and sheared and twisted flux bundles. 

Regarding the magnetic energy, it has been found t hat an active region 



contains enough free magnetic energy to trig ger flares (e.g.. iRegnier et aill2002l : 
Blevbel et al.l 20021 : iRegnier Sz Priestll2007al ). By studying the magnetic energy 



budget before and after a flare, it is difficult to conclude as often the magnetic 
energy released during the flare is in competition with the continuous injec- 
tion of energy from the convection z one as well as the r edistribution of the 
energy inside the volume considered ( Blevbel et al.l l2002l : IRegnier fc Canfield 

I). 



20061 : iThalmann fc Wiegelmannlliooi : ISu et al.ll200 



The magnetic helicity of the magnetic fleld is a quantity more difficult to 
tackle as the knowledge of the vector potential is required inside the coronal 
volume. It has been shown that the magnetic helicity is not a conserved quantity 
in a volume above an active region as magnetic helicity is injected from the 
convection zone and ejected away from the corona. 

In terms of the structure of the corona, the nonlinear force-free field based 
on vector magnetograms has revealed the existence of weakly and highly twisted 
flux bundles in active regions describing solar features such as: 

- Filaments: fllaments are magnetic structures containing cool and dense 
material compared to the coronal environment. In nonlinear force- free 

models, filaments are often identified as we akly twisted fiux bu i idles with 

magnetic di ps in whi ch the plas ma is stored ( Aulanier et al.'l999':'Yan &: Sakurai 
2000: Reg nier Ama ri 2004; Wieeelmann e t al. 2005; Dudik et al. 2001"^ 
Yeates &: Mackavl |2009| ) . Nevertheless observations have shown that ac- 



tive region filaments can be highly twisted as they are subject to kink 
instability; 

Sigmoids: sigmoids are observed in soft X-ray as S or inverse-S shaped 
structures of hot plasma. These structures have been often identified as 



weakly or highly twisted fiux bundles with no magnetic dips (IRegnier &: Amari 
20041 : ICanou et al."2009"; 'Su et al."2009'; 'S avcheva k van Ballegooiienll2009l )" 



As shown by Regnier & Priest (2007b), highly twisted fiux tubes are re- 
quired to store magnetic energy high in the corona; 

Others: other twisted flux bundles are present in magnetic configurations 
with a different amount of twist and/or a different handedness (as both 
signs of currents are observed in a po larity) but which cannot be identified 
to observed features (jRegnier fc Am ari 2004). 



Coronal loops in the core of a ctive regions (observed i n soft X-rays) can carry 
a significant amount of current (jRegnier fc Amari 2004), whilst large loops on 
the edge of active regions (obser ved in EUV at 1-1.5 MK for instance) are close 
to potential magnetic field lines (jPeRosa et al.l [20091 ). 
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It is important to notice that the structure of the magnetic field strongly 
depends on the magnetic field model and on the nature of the active region, 
especially the total magnetic flux a nd the distributio n of polarities on t he phq - 
tosphere (jRegnier fc PriestI l2007bl ^. In particular, H egnier fc PriestI (|2007bl ) 



demonstrated that, statistically, the magnetic field lines are higher and longer 
in a nonlinear force-free configuration than in a potential field one. 

3.3. Full Sun 

The nonlinear force-free description of the whole corona is derived from the flux 
transport model. Starting from a potential field equilibrium, the photospheric 
magnetic distributions are evolved to match the observed synoptic maps and 
the 3D coronal field is thus given by a series of nonlinear force-free equilibria. 
Note that, compare to previous models, the flux transport model does not reset 
the coronal field to a potential field at each time step. As mentioned in the 
previous section, statistically speaking, the field lines in a nonlinear force-free 
model are longer and higher than in a potential field. The consequence is that 
the open magnetic fiux contributing to the fast solar wind is larger in nonlinear 
force-free models than from potential models. The amount of open fiux from 
the fiux transport model is est imated to be one ord er of magnitude larger than 
for the potential field model ( Mackav. D. H. 2010l ). The potential models are 



useful for a qualitative description of the high corona but improved models are 
required to have a better quantitativ e description of the corona. 



As shown bv ICook et al. ( 20091 ). the complexity of coronal magnetic field 



is very low: only few null points (in average, 14 null points) are present in the 
whole 3D corona up to 2.5 solar radii. The time variation of the number of null 
points follows the magnetic cycle. 



4. Discussion 

The nonlinear force- free modelling of the solar corona has becaine a very attrac- 
tive doma in of rese arch ("S chrijver et al.ll2006l : iMetcalf et alll2008l : ISchriiver etaf 



[2008; DeRosa et al. 2009. ) . This physical assumption corresponds to an impor- 
tant step in our understanding of the 3D structure of the solar corona. It corre- 
sponds currently to the state-of-the-art numerical techniques relying on magnetic 
observations. Nevertheless the force- free assumption is debatable especially at a 
time when the new space missions such as the Hinode satellite have significantly 
improved the spatial and time resolutions, and thus show that plasma fiows play 
an important role in the nature of photospheric and chromospheric plasmas. The 
next step to improve this type of modelling based on observations is to consider 
the plasma parameter: the magnetohydrostatic model is a step forward to be 
implemented for future solar missions. As mentioned already, some tentatives 
to model the solar coron a as a magnetohydrostatic equilibrium have attempted 
(|Wiegelmann et all [20071 ). 



The nonlinear force-free models are constrained by the photospheric or 
chromospheric magnetic field. With the large amount of data from Hinode 
or Solar Dynamics Observatory (SDO), it is suggested that more constraints 
should be taken into account to retrieved a more realistic description of the 
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coro nal field. Attempts have been made in se veral papers mentioned above 
(e.g., Ivan BallegooijenI [20041 : 



DeRosa et aLll2009l ^. 
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